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Comment on "Earthquakes Descaled" 

Lindman ef a/. [l|] have used a nonhomogeneous Pois- 
son process with a modified Omori rate, r(t) = dn/dt = 
j*m(1 +t/c)~ p , to model earthquake occurrence. We are going 
to show that contrary to claims in Ref. OL this extremely sim- 
ple description is incomplete in order to explain Bak et a/.'s 
unified scaling law |2fl. 

We generalize Lindman et a/.'s model by introducing an 
r— dependent waiting-time probability density of the form 
D(r\r) °c r^x^ 1 e~ rx l a , which includes the nonhomogeneous 
Poisson process of Ref. QJ], given by y = 1 and a = 1 (both 
parameters linked by normalization). The probability density 
of the waiting times in the Omori sequence, independent of r, 
is given by the mixing of all D(t\r) |3], 


1 fM 

D{x\r m ) = - rD{x\r)p(r)dr 

M J I'm 


(1) 


where p(r) is the density of rates, p(r) °< \dr/dt\~ l = 
C/r l+l l p \ pL is the mean rate of the sequence, fi = / rpdr; 
rM is the maximum rate, corresponding to t = 0; and r m is 
the minimum rate, related to the background seismicity level. 
Note that we have emphasized the dependence on r m . 

Easy to deal with but illuminating is the case y = l/p, 
which yields 


D(l\r m ) °< — . 


(2) 


where the minimum rate r m determines the exponential tail 
of D(z\r m ) for large T, preceded by a decreasing power law 
with exponent 2 — l/p if tm ^> r m . For p = 1 this is in agree- 
ment with the simulation and numerics in Fig. 1 of Ref. Qj]; 
however, it can be shown that the exponent 2 — l/p holds even 
when 7 7^ l/p, which is in disagreement with Lindman et al.'s 
claim of a l/z p decay for p < 1 and 1 /W? for p > 1 . 

Nevertheless, this description totally ignores the spatial de- 
grees of freedom, fundamental in Bak et a/.'s approach. In 
fact, their approach performs a mixing of waiting times com- 
ing from different spatial areas (or cells), which are charac- 
terized by disparate seismic rates. In particular, each area 
will have a different r m , depending on its background seis- 
micity level. This spatial heterogeneity of seismicity can be 
described by a power-law distribution of mean rates R, being 
R the total number of events divided by the total time for a 
given area, see Fig. 1 and Ref. if we assume that the 
minimum rate r m is directly related to the mean rate of the se- 
quence fi, which in turn is in correspondence with the mean 
rate in the area, R, then, p(r m ) °= l/r)~ a and therefore the 
waiting-time probability density comes from the mixing, 


D(t). 


I mil 


r m D(z\r m )p{r m )dr„ 


(3) 


account that C/fl depends on r m , leads to 

D(T) oc 1/t 2+ " for r mm T < 1 < r mM T, (4) 
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FIG. 1 : (color online) Distribution of mean rates R for earthquakes 
with magnitude M > 2 in Southern-California, dividing the area 
(123°W,113°W) x (30°iV,40 o A0 in cells of size L, ranging from 
L = 0.078° to L = 2.5°, and averaging the periods 1984-1992 and 
1993-2001. The distributions are rescaled by L df with df = 1.6. For 
small R/L d f , the data are in agreement with a density p(R) 1 /R ^, 
obtained from a fit. 


which is in disagreement with Lindman et al.'s analysis. 

In fact, the power law for long times [Eq. (0}], was estab- 
lished in Ref. J5| for Southern California, but without relat- 
ing it to the spatial heterogeneity of seismicity. The universal 
value of the exponent 2 + a, found in Ref. J4J] analyzing di- 
verse seismic catalogs, would imply the universality of seis- 
micity spatial heterogeneities. In consequence, Bak et a/.'s 
unified scaling law provides a way to measure these proper- 
ties and is far from being as trivial as suggested by Lindman 
et a/.'s approach. 

The model presented here is still too simple for real seis- 
micity, but provides a clear visualization of its complexity and 
the fundamentals of the unified scaling law of earthquakes. 
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where r m varies between r mm and r,„M- Integration, taking into 
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